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Abstract 

A microscopic model for describing ferroelectric nanoparticles is proposed which allows us to 
calculate the polarization as a function of an external electric field, the temperature, the defect 
concentration and the particle size. The interaction of the constituents of the material, arranged 
in layers, depends on both the coupling strength at the surface and that of defect shells in addition 
to the bulk values. The analysis is based on an Ising model in a transverse field, modified in such a 
manner to study the influence of size and doping effects on the hysteresis loop of the nanoparticles. 
Using a Green's function technique in real space we find the coercive field, the remanent polarization 
and the critical temperature which differ significantly from the bulk behavior. Depending on the 
varying coupling strength due to the kind of doping ions and the surface configuration, the coercive 
field and the remanent polarization can either increase or decrease in comparison to the bulk 
behavior. The theoretical results are compared with a variety of different experimental data. 
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I. INTRODUCTION 



The interesting physical properties and the abihty to store information via a switchable 
polarization have moved ferroelectric materials into the focus of research. One important 
application, the non-volatile memories with ferroelectric capacitor materials are known as 
ferroelectric random access memories (FRAM). The processing issues involved in the high 
density integration process are highly dependent on the ferroelectric and electrode-barrier 
materials. Hence, the selection of materials is a decisive factor in determining the per- 
formance of the device j^. In view of fundamental ferroelectric properties, there are two 
potential ferroelectric materials for FRAM applications, namely Pb(Zr,Ti)Ti03 (PZT) and 
SrBi2Ta209 (SET) P]. They possess a high remanent polarization cx^, low coercive electric 
field Ec, which characterizes the polarization reversal, and low dielectric loss. The use of 
ferroelectric thin films and small particles in high density non- volatile random access mem- 
ories is based on the ability of ferroelectrics being positioned in two opposite polarization 
states by an external electric field j^. It is therefore of great interest to study the thickness 
dependence of Ec for small ferroelectric (FE) particles. This coercive electric field usually in- 
creases significantly with decreasing film thickness or particle size . The strength 
of the coercive field is related to the ease of domain nucleation and domain-wall motion, 
while the permittivity is related to the density of domain walls and their mobility at low 
fields. A diversity of different explanations have been proposed in the past for this size effect, 
including the surface pinning of domain walls and the influence of an internal electric field 
on the domain nucleation in depleted films. Current observations state that the thickness 
dependence of the coercive field in PZT films with metallic electrodes is caused mainly by 
the presence of a non-ferroelectric layer at the film/electrode interface. A parabolic-like re- 
lationship between the coercive field and grain size of ferroelectric PZT films is obtained by 
Liu et al. Nevertheless, the curvature of the parabolic curves varies with the thickness. 
It is demonstrated that the curvature increases as the thickness of the film decreases. The 
influence of thickness on the curvature was explained by the dielectric dependence on the 
thickness. The coercive field of a FE particle is stress sensitive. Lebedev and Akedo ^ 
have shown that the increase of internal compressive stress for thinner PZT films led to the 
increase of the coercive field and breakdown electrical strength, while tensile stress decreases 
ar and Ec. 
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In order to obtain high remanent polarization and low coercive field there are many ex- 
periments with doped ferroelectric thin films and small particles. By adding oxide group 
softeners, hardeners and stabilizers one can modify these materials. Softeners (donors) re- 
duce the coercive field strength and the elastic modulus and increase the permittivity, the 
dielectric constant and the mechanical losses. Doping of hardeners (acceptors) gives hi ghe r 
conductivity, reduces the dielectric constant and increases the mechanical quality factor jlOf . 
Multiple ion occupation of A and/or B sites in ABO3 compounds can affect the lattice 
parameters and the tetragonal distortion (c/a). As a consequence a change of the physical 
properties as the polarization and the phase transition temperature is expected. Direct 
evidence of A-site-deficient SBT and its enhanced ferroelectric characteristics is given by 
Nog uchi et al. The increase of Ba concentration in PLZT ceramics done by Ramam 

Iq and the substitution of La in PZT nanopowders and thin films lead to a marked decrease 



in Tc, ar and Ec 



12, 
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141. An enhancement of dielectric constant and lower and Er 

n 

were observed by addition of PT particles to a PZT matrix |16| . 

There is only a small number of theoretical studies of the dependence of the coercive field 
in ferroelectric small particles. Lo [17j has investigated theoretically the thickness effect in 
ferroelectric thin films using Landau-Khalatnikov theory. Ferroelectric properties such as 
the hysteresis loop, and its associated coercive field and the remanent polarization of various 
film thickness have been numerically simulated. Recently FE thin films has been studied 
in dependence on the film thickness, where the phenomenological Landau theory is applied 
[isj ]. The authors stress that the thickness-dependent coercive field represents the intrinsic 
coercive field. Moreover, they point out to the role of defects observed in real ferroelectric 
thin films. In general the number of defects increases with the film thickness. 
The aim of the present paper is to study the size dependence and doping effects of ferroelec- 
tric small particles using a more refined model, in which the microscopic interaction between 
the constituent parts of the material are explicitly taken into account. A adequate candidate 
is the Ising model in a transverse field (TIM), where ferroelectric properties are described 

nn 

by pseudo-spin operators. The model is quite successful for bulk material [121, 120]. Here 
we extend the applicability of the model to nanoparticles. We utilize a Green's function 
technique to investigate the interacting system, which has to be formulated in real space 
because of the lack of translational invariance. 
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II. THE MODEL AND THE GREEN'S FUNCTION 

Following Blinc and Zeks j3| the FE constituents are described by pseudo-spin operators. 
We consider a spherical particle characterized by fixing the origin at a certain spin in the 
center of the particle. All other spins within the particle are ordered in shells, which are 
numbered by n = 0, 1, A^. Here n = denotes the central spin and n = N represents the 
surface of the system, see Fig. [TJ 




FIG. 1: Array of the ferroelectric nanoparticles composed of different shells. Each sphere represents 
a spin situated in the center, where (a) consists of one central spin plus N = 1 shell, 
(b) A = 2, (c)A = 3, (d) A = 4. 

The surface effects are included by different coupling parameters within the surface shell 
and within the bulk. 

The Hamiltonian of the TIM includes both, bulk and surface properties. It reads 

h = -Iy. '^^^^^^^ - E ^^^^ -f^^Y. (1) 

ij i i 

where and are components of spin-| operators and the sums are performed over 
the internal and surface lattice points, respectively. E is an external electrical field. The 
quantity Jij represents the strength of the interaction between spins at nearest neighbor 
sites i and j. The interaction of spins at the surface shell is denoted as Jij = Jg, whereas 



4 



otherwise the interaction is Jb- In the same manner Vti, and VLg represent transverse fields 
in the bulk and surface shell. It is important to mention that the exchange interaction 
Jij = Jiji — Tj) depends on the distance between the spins, consequently the interaction 
strength is determined by the lattice parameters, the lattice symmetry and the number of 
next nearest neighbors. Starting from Eq. (^, compare one finds immediately that the 
ordered phase is characterized by < S""^ > 7^ and < S"^ > 7^ 0. Therefore it is appropriate 
to introduce a new coordinate system by rotating the original one by the angle 9 in the 
X — z plane. The rotation angle is determined by the requirement < S^' > = in the new 
coordinate system. 

Since the lack of translational invariance the Green's function has to be investigated in the 
real space. The retarded Green's function is defined by 

G',„(t)=«S+(t);S;,(0)» . (2) 

The operators , are the Pauli operators in the rotated system. The equation of motion 
of the Green's function in random phase approximation (RPA) reads 

^Girn = 2{Sl)6im + 21]; siu 9i + flE COS 9i + Jij COS 9i COS 9j (Sj) 

j 

+ J;,sin^;sin0,((5+57) + (5r57))]G;^ 

j 

- \ Y1 '^i^ [^i^ ^^i^ + 2 cos 9j cos 9i{S^Sr)] G.m- (3) 

j 

The poles of the Green's function give the transverse excitation energies. Within the applied 
RPA the transverse spin-wave energy is found as 

e„ = 2Qn sin + ^ Jnj cos 9n cos 9j {Sj)+^E cos 9n , (4) 

j 

where A^' is the number of sites in any of the shells. To complete the soft-mode energy e„ of 
the n-th shell, one needs the rotation angle 6'„, which follows from the condition < S^' > = 0. 
From here the angle is determined by the equation 

— n„ cos 9n + -(^nJu COS 9n siu 9„ + 1^1^^ ^iu ^„ = 0. (5) 

Using the standard procedure for Green's function we get the relative polarization of the 
ra-th shell in the form 

= {S^) = ^ tanh 1^ . (6) 
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III. NUMERICAL RESULTS AND DISCUSSION 



In this section we present the numerical results of the relevant quantities based on Eqs. (jj 
El El)- In particular, we analyze the behavior of spherical FE particles, which are depicted 
in Fig. n In case of BTO the bulk interaction parameter is = 150 K and the transverse 
field is Qb = 10 K. Due to the different numbers of next nearest neighbors on the surface 
the interaction strength J can take different values for the surface, denoted by Jg, and 
for the bulk, labeled as Jb, compare also j21j]. Firstly, let us consider the hysteresis loop 
for different surface configurations represented by the interaction constant Jg at a fixed 
temperature T = 300 K and fixed A^. The results for a particle with = 8 shells are shown 
in Fig. H 

One observes that the coercive field Ec is sensitive to variations of the interaction parameter 
J. If the coupling at the surface is smaller Jg = 50K as that in the bulk Jb = 150 K (dashed 
line), the coercive field Ec and the remanent polarization cr,. are reduced in comparison to the 
case for Js = Jb (solid line). With other words, the coercive field is lowered when the critical 
temperature of the system is decreased due to the smaller Jg value. In the opposite case 
Js > Jb (dotted line) both, the coercive field E^ and the remanent polarization cr^ increase 
compared to Jg = Jb- Because of the enlarged value of Jg the phase transition temperature 
Tc of the small particle will be enhanced. The realization Jg < Jb may explain the decrease 



of the polarization a and the phase transition temperature observed in small BTO- 



22 



and PTO- particles The opposite case, namely a stron ger coupling at the surface, is in 



24| . where the polarization and 



agreement with observations made in small KDP- particles 
the critical temperature increase compared to the bulk material. 

In order to investigate the properties of BTO-type FE small particles we examine in the 
following the case Jg < Jb in more detail. The temperature dependence of the hysteresis 
loop for eight layers is shown in Fig. El 

With increasing temperature the hysteresis loop is more compact and lower, the coercive 
field decreases and even for T > Tc the hysteresis loop vanishes (dashed line). Moreover, we 
have investigated the dependence of the remanent polarization cr^ and the coercive field E^ 
on the particle size A^. The results are depicted in Fig. |31 

Both the field Ec (square) and the polarization cr^ (circle) increase with increasing A^, which 
is in agreement with theoretically finding [18] and the experimental data from ultrathin FE 
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FIG. 2: Influence of the surface coupling strength Js on the hysteresis at fixed temperature T = 
300 K, Jfe = 150 K and iV = 8: Js = 150 K (solid curve), Js = 350K > Jb (dotted curve), 
Js = 50K < Jb (dashed curve); the insert offer the low field behavior . 

BTO- films We found in our approach no rising coercive field with decreasing particle 
size as reported in 

Because there are several experimental indications for a significant influence of doping effects 
on the hysteresis loop, now our model is modified in such a manner to include defects. 
Physically, such defects can be originated by localized vacancies or impurities, doping ions 
with smaller radii and larger distances in-between in comparison to the host material, see 

n 

also Microscopically, the substitution of defects into the material leads to a change of 
the coupling parameter. Within our model let us assume that one or more of the shells are 
composed of defects. The interaction strength between neighbors Jd is altered and in general 
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FIG. 3: Temperature dependence of the hysteresis with N = 8, Jb = 150 K, Jg = 50K: T = 100 K 
(solid curve), T = 300 K (dotted curve ), T = 500 K (dashed curve) . 

different from the surface value Js as well as the bulk one Jb- Furthermore, the defect can 
be situated at different shells within the nanoparticle. In Fig. El we show the influence of 
defects with variable strength Jd on the hysteresis loop. 

A single isolated defect layer offers only a weak influence on the hysteresis curve. Due to 
that we have studied the case that the first five layers of the nanoparticle are defect ones 
(compare also Fig. Q). Here the number of ferroelectric constituents is large enough to 
give a significant contribution to the polarization and consequently to the hysteresis loop. 
Apparently one observes a change of the shape of the hysteresis loop due to defects. In real 
materials one expect Jd < Jb. In that case the coercive field and the remanent polarization 
are reduced (dashed curve) compared to the case without defects Jd = Jb (solid curve). This 
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FIG. 4: Remanent polarization ar (circles) and coercive field Ec (squares) for different layers N at 
temperature T = 300 K and couplings Jb = 150 K, Js = 50K . 

may explain the experimentally observed decrease of T^, Ec and cr^ in small FE particles 
by the substitution of doping ions. The situation is for instance realized by substituting La 



in PTO 12| and PZT [13|, |l^ nanopowders. Although in the majority of real materials 
should be smaller than Jb, we study also the opposite case Jd > Jb- The results are shown 
as the dotted curve in Fig. El The coercive field, the remanent polarization and the critical 
temperature are larger than those without defects (solid curve in (Fig.EJ). Such a situation, 
namely that the defect coupling is stronger than the bulk coupling will be realized, when 
the impurities have a larger radius compared to the constituent ions. The corresponding 
quantities are enhanced in comparison to the bulk value, which is clearly originated by an 
enhanced J^-coupling, which is in accordance with the experimentally observed increase of 
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FIG. 5: Influence of defects on the hysteresis for T = 300 K, Jb = 150 K, Js = 50K . From the total 
number of = 8 layers the first five layers are defect layers: Jd = Jb (solid curve); Jd = 400 K 
(dotted curve); Jrf = 25K (dashed curve). 

Tc, Ec and ar by the substitution of doping ions, such as Bi in SBT [ll| or by increasing 
the Ba contents in PLZT ceramics 3]- This behavior is in contrast to the case Jd < Jb- 
Let us stress that the variation of the interaction strength J can be also interpreted as the 
appearance of local stress, originated by the inclusion of different kinds of defects. The 
case Jd > Jb corresponds to a compressive stress, leading to an enhancement of Ec, which 
has been observed in thin PZT films [9]. The opposite case, i.e. tensile stress, yields to a 
decrease of Ec- 

As visible in Figs. the coercive field and the remanent polarization of the FE particle 

are decreased or increased due to the different interaction strength within the defect shell. 
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Obviously the mentioned quantities should depend on the number of the inner defect shells, 
i. e. on the concentration of the defects. This dependence is shown in Fig. IHl for a particle 
with eight shells. 




2 3 4 5 6 

Number of Defect Layers 



FIG. 6: Coercive field Ec in dependence on the number of defect layers for N = 8, Jt = 
150 K, Js = 50K,T = 300 K and different strength of J^: Jd = 25 K, i. e. Jt, > Jd (squares) and 
Jd = 400 K, i. e. Jd > Jb (diamonds) . 

Notice that for instance Ud = 5 means, that all layers until the fifth layers are defect layers, 
a situation which is precisely considered in Fig. El The squares in Fig. [Ul demonstrates, that 
the coercive field strength decreases with increasing number of defect shells, where we 
assume Jd = 25 K, i.e. Jh > Jd- This result is in reasonable accordance to the experimental 



data given 



= zo IV, I.e. Jh - 
in |liQ,Q. 



A similar result is also obtained for the remanent polarization 



Gt- An increase of the La content in PTO and PZT ceramics decreases the coercive field E^. 
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The opposite behavior is offered as the diamonds in Fig. with Jd = 400 K, i.e. Jd > -h- 
With increasing number of defect layers the coercive field Ec (respectively a^) increases. 
This finding is in a quite good agreement with the experimental data offered by Noguchi et 
al. 11 1 and Ramam et al. jl5|. 

Fig. [7| shows the dependence of the Curie temperature on the number of defect layers 
for a spherical particle. 



800 
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FIG. 7: Dependence of the critical temperature Tc on the number of defect layers Ud for N = 
8, Jfc = 150 K, Js = 50 K and different Jd couplings: Jd = 25 K, i. e. Jf, > Jd (squares) and 
Jd = 400 K, i. e. J5 < Jd (diamonds) . 

A defect coupling smaller than the bulk and surface couplings, Jd < -h (squares), yields 
a decreasing critical temperature Tc with increasing number of defect shells, whereas in 
the opposite case (diamonds) the critical temperature increases with the number of defect 
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layers. The different behavior is controlled by the interaction strength. These results are in 
a qualitative agreement with the experimental data obtained by Noguchi et al. They 
have studied SBT films {T^ = 295 ^C) where the rare earth cations have been substituted 
by La, Ce, Pr, Nd and Sm as well as Bi at the A site (Sr site) with Sr vacancies. The La 
modification induces a soft behavior (lowering of and Pr), while a large amount of Nd and 
Sm substitution results in a very high Ec (hard). So the modified behavior is the result of 
defect engineering of both Sr and oxide vacancies. Likewise other bismuth layer-structured 
FE (BLSF) reveals that is strongly influenced by the radius rj of A-site cations. BLSF 
with smaller A-site cations (Ca^^) tends to a higher Tc (420 °C). The inclusion of the same 
amount of Ba^"*" (with larger radius) leads to a relaxation of FE distortions and results in a 
decrease of to 120*^0. The substitution of La with concentration x = 0.5 gives rise to a 
marked decrease in Tc to 180 °C, because the induced A-site vacancies weaken the coupling 



between neighboring BOg octahedral 



261 El- 

This result corresponds in our approach 



when the interaction parameter Jd is reduced, i. e. Jd < Jb- In La-modified PbTiOa the 
ransition temperature Tc decreases significantly, too, with an increase of the La content 
28l |. In Bi-SBT the critical temperature Tc increases strongly to 405^0 (x=0.2) Q- The 
increase in Tc by Bi substitution tends to the opposite behavior as that observed in La- 
SBT. This is governed by the bonding characteristics with oxide ions. The influence of the 
orbital hybridization on Tc is very large, and Bi substitution results in a higher Tc Q- This 
experimental observation is in a sufficient good agreement with our analytical finding when 
the interaction strength satisfies Jd > Jb- 



IV. CONCLUSION 



Because ferroelectricity at a nanoscale has emerged as a fertile ground for new physical 
phenomena we have analyzed in this paper FE nanoparticles based on a microscopic model. 
To this aim the Ising model in a transverse field is modified in such a manner that the 
effects of surfaces and defects are taken into account. Eventually, these properties are 
manifested within our model by introducing different microscopic coupling parameters for 
the constituents of the FE material. Due to the broken translational invariance the Green's 
function technique has to be formulated in real space. As the result of the equation of motion 
method for the Green's function we find the polarization of the nanoparticle depending on 
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an external electric field, the temperature, the defect position and concentration as well as 
the size of the particle. In particular, we are interested in the coercive field Ec, which is very 
sensitive to the interaction strength at the surface Jg and the presence of defects manifested 
by a microscopic coupling J^- In case of a defect coupling strength that is lower than the bulk 
one, Jd < Jb, the coercive field, the remanent polarization Cr and the Curie temperature 
Tc are reduced in comparison to the case without defects. This theoretical finding gives 
an explanation of the experimentally observations in FE small particles by substituting of 
doping ions, such as La in PTO [11], PZT [12,13] nanopowders. Contrary we get in the 
opposite case with Jd > Jb (for example when the impurities have a larger radius compared 
with the constituent ions), that the coercive field, the remanent polarization and the critical 
temperature are larger than the case without defects. This realization is in accordance 
with the experimental data showing an increase of Tc, Ec and ar by the substitution of 
doping ions, such as Bi in SBT [10] or by increasing Ba contents in PLZT [14] ceramics. 
The dependence on the particle size is also discussed. Our results give strong indications 
that microscopic details of the interaction within the FE nanoparticles are essential for the 
macroscopic behavior of such a quantity as the hysteresis loop, whose shape seems to be 
important for practical applications of FE. In this paper we have demonstrated that one of 
the standard model for describing ferroelectric systems, namely the modified Ising model in 
a transverse field, is able to give explanation for experimental observations. 
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